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Abstract 

In the present work, the Brans-Dicke (BD) theory of gravity is taken as a possible theory of 
k-essence. Then starting with the (already known) Brans-Dicke-Schwarzschild solution which can 
represent the gravitationally bound static configurations of the BD scalar k-essence, issues like 
whether these configurations can reproduce the observed properties of galactic dark matter haloes 
have been addressed. It has been realized that indeed the BD scalar k-essence can cluster into dark 
matter halo- like objects with flattened rotation curves while exhibiting a dark energy- like negative 
pressure on larger scales. 
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I. Introduction 



Currently, perhaps the most fashionable candidates for the unified model of dark matter 
and dark energy jjj with non-trivial dynamics are quintessence and k-essence. The main 
difference between the two models is that the quintessence models |2| involve canonical 

n 

kinetic terms and the sound speed of (r s — 1 while the k-essence models j2| employ rather 
exotic scalar fields with non-canonical (non-linear) kinetic terms which typically lead to 
the negative pressure. And the most remarkable property of these k-essence models is 
that the typical k-essence field can overtake the matter energy density and induce cosmic 
acceleration only at the onset of the matter-dominated era and particularly at about the 
present epoch. These models are also expected to provide a successful explanation of the 
phenomena associated with the dark matter. In the present work, we take the Brans-Dicke 
(BD) theory of gravity |4[ as a possible k-essence theory since it involves probably the 
simplest form of such non-linear kinetic term for the (BD) scalar field. Besides, the BD 
scalar field (and the BD theory itself) is not of quantum origin. Rather it is classical in 
nature and hence can be expected to serve as a very relevant candidate to play some role in 
the late-time evolution of the universe such as the present epoch. Indeed, the BD theory is 
the most studied and hence the best-known of all the alternative theories of classical gravity 
to Einstein's general relativity lal- This theory can be thought of as a minimal extension of 

n n 

general relativity designed toproperly accomodate both Mach's principle p, |6| and Dirac's 
large number hypothesis j^, H|. Namely, the theory employs the viewpoint in which the 
Newton's constant G is allowed to vary with space and time and can be written in terms 
of a scalar ("BD scalar") field as G = 1/$. As a scalar-tensor theory of gravity, it involves 
an adjustable but undetermined "BD-parameter" to and as is well-known, the larger the 
value of u, the more dominant the tensor (curvature) degree and the smaller the value of 
u), the larger the effect of the BD scalar. And as long as we select sufficiently large value 
of u, the predictions of the theory agree perfectly with all the observations/experiments to 
date [5]. For this reason, the BD theory has remained a viable theory of classical gravity. 
However, no particularly overriding reason thus far has ever emerged to take it seriously 
over the general relativity. As shall be presented shortly in this work, here we emphasize 
that it is the existence of dark matter (and dark energy as well, see |7[) that puts the BD 
theory over the general relativity as a more relevant theory of classical gravity consistent 
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with observations that have so far been unexplained within the context of general relativity. 



II. Haloes of BD scalar k-essence 



In general, the Brans-Dicke theory of gravity is described, in the absence of ordinary 
matter, by the action 



S 
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where $ is the BD scalar field representing the inverse of Newton's constant which is allowed 
to vary with space and time and u is the generic dimensionless parameter of the theory. 
Extremizing this action then with respect to the metric g^ v and the BD scalar field $ yields 
the classical field equations given respectively by 
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Note that here in the present work, we are interested in the role played by the BD scalar 
field (i.e., a k-essence) as a dark matter particularly in forming galactic dark matter haloes 
inside of which the well-known rotation curves have been observed. Since the galactic dark 
matter haloes are roughly static and spherically-symmetric, we first should look for such 
dark matter halo-like solution of these BD field equations. Interestingly enough, the Brans- 
Dicke-Schwarzschild (BDS) spacetime solution to these vacuum BD field equations that 
happens to meet our above-mentioned needs has been found some time ago |S| in rather a 
theoretical attempt to construct non-trivial black hole spacetime solutions in BD theory. To 
summarize, the BDS spacetime solution that can be obtained by setting a = e = in the 
Brans-Dicke-Kerr- Newman (BDKN) solution in eq.(ll) of |8( takes the form 
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^2/(2^+3) sin 4/(2a,+3) g 



where A = r(r — 2M). A remarkable feature of this BDS solution is the fact that, unlike 
the Schwarzschild solution in general relativity, the spacetime it describes is static (i.e., 



non-rotating) but not spherically-symmetric. Of great interest in this earlier construction 
was the realization that non-trivial black hole solutions different from general relativistic 
solutions could occur in this BD theory for the generic BD-parameter values in the range 
—5/2 < uj < —3/2 ^]. In the present study, however, since we are interested in the 
galactic halo-like configuration, we do not want this BDS solution to become "black" and 
this amounts to considering the BDS solution having the value of cj-parameter well outside 
this range. Besides, we are only interested in whether the self-gravitating k-essence, i.e., 
the BD scalar field can generally cluster into dark matter halo-like objects which would be 
the gravitationally bound static solution configurations of super-galactic scale (i.e., the large 
but fmite-r behavior). Therefore, the peculiar microscopic geometrical nature of this BDS 
solution such as the issue of regularity of the potential Killing horizon (i.e., the fmiteness 
of the invariant curvature polynomials there) addressed in |8j or that of seemingly failure of 
asymptotic flatness and internal infinity nature of the symmetry axis discussed in P| are all 
irrelevant for the present purposes. 

Therefore first, it appears that the BD scalar k-essence can indeed cluster into halo-like 
configurations as it can be represented by the BDS solution. Our natural next mission 
is then to ask whether these configurations really can reproduce the properties of dark 
matter haloes, namely if our BD scalar k-essence model for dark matter can reproduce the 
flattening of the rotation velocity curves inside these halo configurations consistent with the 
observations. Thus we now attempt to obtain the rotation curves in our BD scalar k-essence 
halo. Since henceforth we need concrete "numbers" , we now restore both Newton's constant 
Go and the speed of light c in order to come from the geometrical unit (Go = c = 1) back 
to the CGS (or MKS) unit. Then the energy-momentum tensor of the BD scalar field given 
earlier in eq.(2) now takes the form 

T ™ = 8^ [^ (V ^ $V ^ " ^V a $V a $) + ^(V,V,$ - ^V Q V Q $)] (4) 

and in the BDS solution in eq.(3) above, we should replace M —>■ GqM/c 2 = M where Go 
denotes the present value of the Newton's constant. Apparently, (4) has the dimension of the 
energy-momentum density in the CGS unit, (erg /cm 3 ). We now turn to the computation 
of energy density profile and (anisotropic) pressure components of the k-essence playing the 
role of the dark matter by treating the BD scalar field as a (dark matter) fluid. The BD 
scalar field fluid, however, would fail to be a "perfect" fluid as can readily be envisaged 
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from the fact that the associated BDS solution configuration is not spherically-symmetric. 
Namely, its pressure cannot be "isotropic", i.e., P r ^ Pe ^ P<p. Such fluid may be called 
imperfect fluid due to the anisotropic pressure components and as such its stress tensor can 
be written as 

/ -c 2 p 
P r 
T e r 
V 

And it is to be contrasted to its counterpart of the usual perfect fluid with isotropic pressure 
given by the well-known form, T*f = P5% + (c 2 p + P)U^U U = diag(-c 2 p, P, P, P) where 
U a = dX a /dr (with r being the proper time) denotes the 4- velocity of the fluid element 
normalized such that U a U a = —1. Note that in addition to the diagonal entries representing 
the (anisotropic) pressure components T\ = Pi (where no sum over i), there are off-diagonal 
entries T e r , Tf; representing a shear stress which also results from the failure of spherical 
symmetry. Thus by substituting the BDS solution given in eq.(3) into the BD energy- 
momentumm tensor in eq.(4) and then setting (4) equal to (5), we can eventually read off 
the energy density and the pressure components of the BD scalar field imperfect fluid to be 
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Note that the off-diagonal components Tg, Tf! are odd functions of 9 while the diagonal 
components (p, P r , Pq, P^) are even functions of the polar angle under 9 — > (7r — 0). As 



a result, the off-diagonal components vanish (i.e., no shear stress survives) if we average 
over this polar angle to get a net stress. Thus, first the equation of state of this BD scalar 
k-essence fluid forming a galactic halo is given by 
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c 2 P ~ 2(u + l){(r-M) 2 + Acot 2 #} - {2u + 3)M(r — M) 

where P = P r . Namely, P = w(r,6)c 2 p with w(r,6) ~ 0(1) meaning that this k-essence 
fluid is essentially a barotropic fluid but with "position-dependent" coefficient w(r,8). Note 
that although the BD scalar k-essence is a candidate for dark matter, it is not quite a 
dust. In principle, the speed of sound in this BD scalar field fluid can also be evaluated 
via c 2 = dP/dp but we shall not discuss it in any more detail in this work. We are now 
ready to compute the behavior of rotation curves in the outer region (i.e., at large but 
finite-r, say, r » GqM/c 2 ) of our BD scalar k-essence halo. To be more precise, for a 
galaxy of typical (total) mass M ~ 10 11 M Q , the outer region of its dark matter halo, say, 
r ~ 10(kpc) ~ 10 23 (cm) is much greater than G M/c 2 ~ 10 16 (cm) by a factor of 10 7 . Thus 
to this end, we first approximate the expressions for the energy density and the (radial) 
pressure of the k-essence given in eq.(6) for large-r. They are 

° 2 8 ( U + l l \ A -2/(2u;+3) j -4/(2uH-3)0 (8) 
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Note that in the above approximations and in the discussions below, it was and it shall be 
assumed that the metric function A = r(r — 2M) ~ r 2 for large-r. It is interesting to note 
that as a "k-essence" constituting a dark matter halo, the energy density p of the BD scalar 
field is almost certainly positive everywhere (i.e., for both small and large-r). In the mean 
time, its (radial) pressure P particularly at larger scale (i.e., for large-r) turns out to be 
negative although its sign appears unclear at small scale (i.e., for small-r). 
Finally, we are ready to determine the rotation curve inside our BD scalar k-essence halo. 
First in the most naive sense, the apparent rotation velocity of an object at radius r from 
the galactic center is given by the Kepler's third law, v 2 = G M(r)/r. Thus for our 
case, using the BD scalar k-essence energy density profile given earlier, we have M(r) = 
C # i: d6 / r dr^g rr g ee g H p{r, 6) = (2c 2 /G ) [(w + l)/(2w + l)(2w + 3)] /(c ,) rr - 2 /( 2 -+3) 
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and hence 



2 . . G M(r) 2 2(u + l) 



r (2a; + l)(2u; + 3) 

where f(u) = g d9 s i n -[i+2/(^+3)] = 2 /^^[l - ^-(a^y^+s) with ^ tf L ( Note 
here that the integration over the polar angle 9 starts not from but from e << 1 as the 
symmetry axis 9 = of the BDS solution in eq.(3) possesses an internal infinity nature, 
namely, the symmetry axis is infinite proper distance away as discussed carefully in |9j.) It 
has been known for some time that in order for the BD theory to remain a viable theory 
of classical gravity passing all the oblational/experimental tests to date, the BD u,- 
parameter has to have a large value, say, \u\ > 500 [5]. In our previous study jSj, in the 
mean time, it has been realized that the static solution to the vacuum BD field equations 
given in eq.(3) above can turn into a black hole spacetime for —5/2 < u < —3/2. Thus 
now for \u\ > 500, the same static solution eq.(3) we are considering represents just a halo- 
like configuration with regular geometry everywhere (i.e., having no horizon) which is static 
but not exactly spherically-symmetric (note that the galactic haloes are also believed to be 
nearly spherically-symmetric but not exactly). Thus if we substitute a large-u; value, say, 
uj ~ 10 6 into eq.(9) above, evidently M(r) ~ r and hence we get 

v(r) ~ 100(Jfem/s) x r ~ (1/1 ° 6) (10) 

since for to ~ 10 6 , f(u) — 0(1). Namely for this large-u; value, the rotation curve gets 
flattened out as r~( 10 6 - ) ~ constant and its magnitude becomes several hundred (km/s). 
Indeed, this is in impressive agreement with the data for rotation curves observed in spi- 
ral/elliptic galaxies with M/L ~ (10— 2O)M /L and in low-surface-brightness (LSB)/dwarf 
galaxies with M/L ~ (200 — 6OO)M /L Q (where M/L denotes the so-called "mass-to-light" 
ratio given in the unit of solar mass-to-luminosity ratio exhibiting the large excess of dark 
matter over the luminous matter) Rotation curves are observed usually via the mea- 
surements of the Doppler shift of the 21cm emission line from neutral hy drog en (HI) for 
distant galaxies and of the light emitted by stars for nearby galaxies (jjj, [l2j. It is also 
interesting to note that this behavior of the rotation curve in our BD theory k-essence dark 
matter halo model is independent of the mass of the host galaxy as it should be. Namely, 
this behavior of the rotation curve comes exclusively from the nature of the dark matter, 
i.e., the BD scalar field k-essence. We also point out that even if we employ more careful 



expression for the rotation velocity curve involving the Doppler shift of light emitted by the 



orbiting objects (assuming that the 
v 2 (r) = G M(r)/r + 4vrr 2 G P/c 2 



^-essence halo is almost spherically-symmetric), namely 
(with P being the radial pressure given in eq.(8) 
above), the conclusions above remain the same. 

Next, the equation of state in eq.(7) of this BD scalar k-essence becomes, in the outer region 
of the galactic dark matter halo (i.e., at large-r), 

2(u - 1) cos 2 6> + sin 2 ' 



W ~ — ; r (11) 

2(w + l) V ; 

which is obviously negative due to the negative pressure (and still positive energy density) 
in this outer region. Moreover, for the large-u; value, i.e., uj ~ 10 6 for which the rotation 
curve gets flattened out that we just have realized, this equation of state at large-r further 
approaches w ~ — cos 2 9 ~ —0(1). (Incidentally, it is interesting to note that in the vicinity 
of the equatorial plane 9 = n/2, w = 0, namely, the BD scalar k-essence behaves like nearly 
a dust.) This observation is particularly interesting as it appears to indicate that the BD 
scalar k-essence we are considering possesses dark energy-like negative pressure on larger 
scales. And this observation is indeed consistent with our previous study |7[ that on the 
cosmological scale, the BD scalar field does exhibit the nature of dark energy possessing the 
negative pressure. 

III. Concluding remarks 

In the present work, starting with the (already known) BDS solution which can represent 
the gravitationally bound static configurations of the BD scalar k-essence, issues like 
whether these configurations can reproduce the observed properties of galactic dark matter 
haloes have been investigated. It has been realized that indeed the BD scalar k-essence can 
cluster into dark matter halo-like objects with flattened rotation curves while exhibiting a 
dark energy-like negative pressure on larger scales. Thus to conclude, from this success of 
"BD scalar field as a k-essence" to account for the asymptotic flattening of galaxy rotation 
curves while forming galactic dark matter haloes plus the original spirit of the BD theory in 
which the BD scalar field is prescribed not to have direct interaction with ordinary matter 
fields (in order not to interfere with the great success of equivalence principle), we suggest 
that the Brans-Dicke theory of gravity is a very promising theory of dark matter. And this 



implies, among others, that dark matter (and dark energy as well, see j^l) might not be 
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some kind of unknown exotic "matter" , but the effect resulting from the space-time varying 
nature of the Newton's constant represented by a (k-essence) scalar field. Even further, this 
successful account of the phenomena associated with the dark matter of the present universe 
via the BD gravity theory might be an indication that the truly relevant theory of classical 
gravity at the present epoch is not general relativity but its simplest extension, the Brans- 
Dicke theory with its generic parameter value u ~ 10 6 fixed by the dark matter observation ! 

This work was financially supported by the BK21 Project of the Korean Government. 
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